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THE LONGITUDINAL STABILITY OF ELASTIC SWEPT WINGS AT SUPERSONIC SPEED 

By C. W. Frick and R. S. Chubb 


SUMMARY 

The longitudinal stability characteristics of elastic swept 
wings of high aspect ratio experiencing bending and torsional 
deformations are calculated for supersonic speed by the appli- 
cation of linearized lifting-surface theory. A parabolic wing 
deflection curve is assumed and the analysis is simplified by 
a number of structural approximations. The method is thereby 
limited in application to wings of high aspect ratio for which 
the root effects are small. Expressions for the lift, pitching- 
moment, and span load distribution characteristics are derived 
■in terms of the elastic properties of the wing; namely, the design 
stress, the modulus of elasticity, the shearing modulus, and the 
maximum design load factor. The analysis applies to wings 
with leading edges swept behind the Mach lines. In all cases, 
however, the trailing edge is sonic or supersonic. Application 
of the method of analysis to wings with leading edges swept 
ahead of the Mach lines is discussed. 

The results of numerical calculations for a wing of aspect 
ratio 3.2 and 60° sweeplack are presented for a Mach number 
of 1M4 an d for incompressible flow. The effects of wing 
elasticity on the lift-curve slope, moment-curve slope, and neu- 
tral-point position are shown. The results indicate that the 
primary variable involved in aeroelastic phenomena is the 
dynamic pressure and that the influence of the flight Mach 
number is small for wings swept behind the Mach lines. 

INTRODUCTION 

la reference 1, R. T. Jones has shown that supersonic 
flight may be attained with a reasonable degree of efficiency 
through the use of swept wings of high aspect ratio. The 
use of sweepback, however, involves many problems of 
stability and control, not the least of which are associated 
with the aerodynamic effects of the elastic deformation of 
the airplane structure. In particular, the longitudinal 
stability of the aircraft may be affected to a large degree 
since the bending and torsional deformations of the wing 
may shift the center of pressure of the lift forward an appre- 
ciable distance. 

These aeroelastic phenomena occur under those flight 
conditions where the magnitude and/or the spanwise varia- 
tion of the elastic deformation of the wing varies with 
angle of attack. Aeroelastic effects may therefore occur 
either in accelerated flight at constant dynamic pressure or, 
under certain conditions, in steady level flight with varying 
dynamic pressure. In the latter case, if the loading due to 
twist or camber is different than the loading due to change 
of angle of attack, the trim change due to elastic deformation 
of the wing in steady level flight varies with the dynamic 
pressure and influences the stability of the airplane as 
indicated by the position of the control stick as a function 
of airspeed. 1 


In solving aeroelastic problems, since the interrelation of 
the structural and aerodynamic characteristics of the wing 
results in mathematical complexity, it is us ually necessary 
to compromise to some extent either the structural or the 
aerodynamic aspects of the problem to obtain a solution. 
In the present analysis, the structural characteristics of the 
wing are compromised to the extent that the form of the 
deflection curve is assumed. This assumption permits the 
application of supersonic lifting-surface theory to the deter- 
mination of the load distribution, the lift, and the pitching- 
moment characteristics of elastic wings. Additional analysis 
is necessary to determine whether it is better to use more 
rigorous aerodynamic theory in aeroelastic computations, as 
in the present report, or to use a more complete structural 
theory as in recent work by Diederich (reference 2) and 
Miles (reference 3). 

SYMBOLS 


X, ,y l Cartesian coordinates 

x,y transformed Cartesian coordinates in terms of the 
semispan dimension, s 

£, 1 ? x,y coordinates of the apex of any superposed lifting 
sector 

s distance in the y x direction from the root section to the 
intersection of the flexural axis and the tip Mach 
cone 

s' distance along the flexural axis from the root- section 
to the intersection of the flexural axis and the tip 
Mach cone 

y' distance measured from the- root section along the 

flexural axis 

7 spanwise distance in y direction from the root section 
to the center of load on the half wing 

S' wing area 

A taper ratio, ratio of tip chord to root chord 

c a , average chord 


c 


mean aerodynamic chord 



c local chord parallel to the plane of symmetry 

c 0 root chord parallel to the plane of symmetry in terms 

of the span dimension, s 
AR aspect ratio 

A angle of sweepback of the flexural axis 
6 slope of the flexural axis in a vertical plane passing 
through the flexural axis 
n maximum load factor 

M bending moment at any point on the flexural axis 
M T bending moment at the root section of the wing beam 


1 This particular aeroelastic characteristic is not considered in the present report which Is 
concerned primarily with accelerated flight. Further, the wing is considered to be weightless 
so that the ameliorating influence of the distributed mass of the wing is not taken into account 
in estimating the aeroelastic characteristics. 
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torsional moment at any point on the flexural axis 
torsional moment at the root section of the wing beam 
modulus of elasticity for the wing beam material 
shearing modulus for the wing beam material 
moment of inertia of the wing beam 
torsional stiffness constant 

distance between the flexural axis and the center of 
pressure of the sectional lift in terms of the local 
chord 

maximum design stress 

maximum thickness of the wing at the root section 
angle of attack of the root section of the wing 
incremental angle of attack at any spanwise station of 
the wing 

angle of attack of the wing section at any spanwise 
station 

angle of attack of the root section at jvhich maximum 
load factor is developed 

■yjM 2 — 1 where M is the free-stream Jvlach number 
j3 times the cotangent of the angle of sweepback of the 
wing leading edge 

/3 times the cotangent of the angle of sweepback of the 
wing trailing edge 

/} times the cotangent of the 'angle of sweepback of a 
ray from the apex of any superposed lifting sector 
complete elliptic integral of the second kind with 
modulus (Vl — m 2 ) 
airplane weight 

wing loading 


dynamic pressure 



where p is the mass density 


and V the velocity of the free stream 


lifting pressure coefficient 

load per unit span 
section lift coefficient 
lift 

lift coefficient 

lift coefficient at maximum load factor 
section pitching moment of a wing section about the 
apex of the wing 

pitching-moment coefficient about the apex of the wing 
in terms of the mean aerodynamic chord and the 
wing area 

the rate of change of lift coefficient with the angle of 
attack of the root section 

the rate of change of pitching-moment coefficient with 
the angle of attack of the root section 
the rate of change of pitching-moment coefficient with 
the lift coefficient 



ANALYSIS 

WING WITH A SUBSONIC LEADING EDGE 


In the following analysis, for convenience, the aerodynamic 
loading due to bending and that due to torsion are first 
treated separately. Expressions for the combined effects ol 
bending and torsion are derived later. 


Bending. — The aerodynamic twist 2 due to bending of a 
stream wise section of an elastic swept wing under accelerated 
flight conditions is a function of the applied load and the 
elastic characteristics of the wing beam. In order to arrive 
at a solution for the aerodynamic properties of the wing 
without becoming involved in laborious graphical analysis, 
some simplifying approximations must be made regarding 
the elastic properties of the wing. 

In a strict sense, a swept wing of conventional structural 
design cannot be considered to have a flexural axis. For 
wings of high aspect ratio, however, it will be assumed that 
a flexural axis exists, since this assumption permits the use 
of simple beam theory and introduces only a small conserva- 
tive error. 

For the purpose of analysis, the root section of the wing 
beam is assumed to be the extension of the wing beam on a 
plane perpendicular to the flexual axis and passing through 
the intersection of the flexural axis and the streamwise root 



section. (See fig. 1.) This simplification of the beam analy- 
sis is similar to that of reference 4. The length of the wing 
beam s' is taken as the distance along the flexural axis from 
the root to the intersection of the flexural axis and the tip 
Mach cone. The semispan s of the wing is taken as extend- 
ing from the root section to the intersection of the flexural 
axis and the tip Mach cone in a direction perpendicular to 
the plane of symmetry. The portion of the wing lying 
within the tip Mach cone is ignored since, as shown in 

a The change in camber of the airfoil sections due to the distortion of the wing surface is, 
of course, ignored. 
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Figcbe 2.— Coordinate system for calculation of characteristics of elastic wings. 


reference 5, very little load is carried in this region and the 
analysis is thereby simplified. 

The coordinate system is selected as shown in figure 2. 
The origin of the coordinate system is placed at the apex of 
the wing, the positive branch of the ari axis lying downstream. 

The mathematical treatment may be made less tedious by 
transforming and nondimensionalizing the coordinates so 
that in the following analysis 



root chord 

c °~ i 

In general, at both subsonic and supersonic speeds, selec- 
tion of the wing plan form for low drag leads to a combination 
of spanwise loading and spanwise distribution of the bending 
resistance in the wing beam such that the wing deflection 
curve is essentially parabolic. (The ratio of M to I is con- 
stant across the span.) The deflection curve deviates 
appreciably from a parabola only if the aeroelastic effects 
experienced by the wing are very large. 


Since the deflection curve of the flexural axis is assumed to 
be parabolic, the slope of the flexural axis is 


„ M , 

d ~EI y 


where y' is measured along the flexural axis. 

The incremental angle of attack of streamwise sections of 
the wing is related to the slope of the flexural axis as 

a T — — 6 sin A 

The slope of the flexural axis in nondimensional transformed 
coordinates may be written as 

.Ms 

6 ~EI$ cos A y 

The incremental angle of attack of any streamwise section of 
the elastic wing is then 

M s . , 

ok? — -ptj g V tan A 

and the total angle of attack of any streamwise section is 

jfi-py tan A (1) 


where a is the angle of attack of the root section of the wing. 
Equation (1) gives the magnitude and distribution of twist 
across the span of the wing if the magnitude of M/EI is 
known. 

The distribution of pressure over the elastic wing due to 
twist may be determined by applying known conical-flow 
solutions for supersonic flow. In the linearized theory, the 
principle of superposition of various solutions may be used to 
satisfy the particular boundary conditions of the problem. 
For the elastic wing, the flow field may be considered to 
consist of the superposition of two distinct flow fields : 

1. The flow about a flat rigid wing at an angle of attack 
equal to the angle of attack of the root section. 

2. The flow about a twisted wing for which the angle of 
attack at the root is zero. 

The solution for the first flow field is given in references 6 
and 7; the second flow field can be obtained by determining 
the solution for a differential twist da, at one station and 
integrating this solution across the span. 

The solution for the pressure distribution corresponding 
to a differential twist must meet the following boundary 
conditions (fig. 3): 

1. Outboard of the station of twist, the angle of attack 
must be constant and equal to the differential twist. 

2. Inboard of the station of twist, the angle of attack of 
the surface must be zero. 

3. Between the swept leading edge and the Mach cone, 
no lifting pressures may exist. 
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Figure 3. — The boundary conditions and pressure distribution for the conical-flow solution for differential twist. 


The conical-flow solution corresponding to these boundary 
conditions is that for a special lifting sector given by Lager- 
strom in reference 8 and is expressed in the notation of the 
present report as 

Ap_8 a m 3/2 j l+t /g\ 


of the wing r)—0 to the value of v—Vo corresponding to the 
last superposed sector, the Mach cone of which encompasses 
the point x,y under consideration. The value of ij 0 is found 
by placing t equal to -1 and solving for ij. 

570 = ^k <*+*> 


where t defines a ray from the apex of the sector. 

Figure 3 shows both a sketch of the boundary conditions 
to be met by this solution and a plot of the pressure distri- 
bution given by equation (2). 

The induced pressure resulting from twist due to bending 
of the elastic wing may be found by integrating across the 
span of the wing. This integration corresponds to the 
superposition of an infinite number of the lifting sectors 
along the span, each sector having an infinitesimal angle of 
attack da s . 

The pressure due to twist is then given by_ 

8 m? n Ms 


/Ap\ 

\ 2 )t 


/SV m+ 1 El 


tan A 


C^Jl+ldr, 

Jo V m—t 


where 


.y —V = m(y—rj) 

x—% mx—y 


The x and y coordinates of the apex of any superposed sector 
are £, ij. 

The integration must be carried out from the root section 


The integration yields at any point x,y the pressure due 
to twist 

/Ap\ 16 m 5/2 t , Ms , , ^ I x+y 

\ 2 )t~ 3/3 V (m+1) 2 fcan A El (a!+ ^V mx-y ( } 

To this expression must be added the conjugate term due to 
the elastic deformation of the opposite wing panel. The 
conjugate term may be obtained by substituting -y for y. 
Then 

/Ap\ _ 16 Ms T ... / x+y 

V g L . (w+i? ei l {x+v) y^q} + 

^-»>ViS+v] (4> 

It should be noted that the addition of the conjugate 
terms adds some very small lifting pressure in the region 
between the wing leading edge and the Mach cone where no 
lifting pressure may exist. These pressures may be canceled 
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by the superposition of constant lift sectors as noted in 
reference 9. Since these extraneous pressures, on the 
average, amount to about 3 percent of the average pressure 
coefficient over the adjacent wing surface and, since elimi- 
nation of those pressures would change the pressures over 
the surface only about one-half of 1 percent, it seems that 
in view of the additional complication involved the can- 
cellation of these pressures is unwarranted. 

The total lifting pressure for the elastic wing at an angle 
of attack is then obtained by adding to equation (4) the 
solution for the flat lifting wing. For the elastic wing, then 

A p 4mV 16 m s/2 x MsC . . . / x-\-y 

IT ~~02 V m 2 — (y jx ) 2 “ 1 

+(x-y)J— w-l (5) 

Examination of this equation shows that the relationship 
between M/EI and a must be established before the pressure 


distribution can be calculated. Since for wings with para- 
bolic deflection curves the maximum stress occurs at the 
point of maximum thickness, usually the root, the maximum 

stress occurring at maximum load factor is 
* 

<k 

<r ““~ \l) n 2 

and, since the bending moment at any point on the span is a 
linear function of the angle of attack, 

H ^*&jnax 

I d r a* 

where c mai is the design stress at maximum load factor, d, 
is the thickness of the root section and a* is the angle of 
attack at maximum load factor; an expression for a n is 
derived later. 

The equation for the pressure distribution may then be 
written as 


Ap 4 mV 

2 _ 02 ■Jm 2 — (yjx ) 2 


32 m 5/2 

30V (m+1) 2 tan A 


E d, 




g— y 

mx+yj 


(6) 


The load per unit span can be obtained from an integration 
of equation (6) with respect to x along any streamwise station 
(y-constant) , 

V+m t c 0 


1 r f, 

2 S X \ 


Ag\ 

2 /ff=const. 


dx 


The integration is carried out from the leading edge of the 
wing, x=— to the trailing edge and yields 


m 

_Z 4 m?sa 
2 


7/lt 


02 


My)- 


32 m 5/2 ffm,, s* a , . . 

'30V (m+l) stanA E d r ccj i(y) (/) 


The functions /i(y) and Mv) 816 given in the appendix since 
they are somewhat unwieldy. 

The lift coefficient may be obtained by an integration of 
equation (7), spanwise from root to tip. 


\ iy 


The integration yields 3 


0 &-§-[- 


2s 2 T 4 mV 
02 


F\- 


32 m s/2 
'30V (m+1) 2 


tan A ■ 


E d. 


r~F 2 1 (8) 

r a n J 


The constants Fi and F 2 are given by equations in the 
appendix. 

This equation may be used to determine the angle of attack 
at maximum load factor a n which is needed in the foregoing 
equations: 


0"2 r<s n 

““ 4m 2 F 1 |_2s 2 L 


32 


m 


5/2 


>tan A- 


|F 2 ] (9) 


30 V (m+1) 2 E dr 

The pitching-moment characteristics of the elastic wing may 


3 It may be noted that the ratio s 2 fS is essentially the same as one-fourth, aspect ratio and 
that the parameter a}d r is directly related to tffdt, a common structural criterion. 


be determined by an integration of the pressure distribution 
given by equation (6) . 

For any spanwise station, the section pitching moment 
about the apex of the wing is 



v+mco 


d - m 


x dx 


This integration yields 



32 m 5/2 

30 V (m+1) 2 


tan A 


v max $ cc. 

E drOn 



The functions My) and My) are given in the appendix. 

The total pitching-moment coefficient about the apex of 
the wing in terms of the mean aerodynamic chord is found 
by integration across the span, 



BC - 2s s [~4toV-, 

P Lm ~ ScL 02 


m* 


5/2 


32 

30 V (m+l) ! 


;tan A- 


E d -. 




or 


BC =— 

Set I 


m 2 F s 32 


IK 


5/2 


02 30V(m+l) : 


;tan A- 


E d. 


<12> 


The constants F z and F t which are functions of the aspect 
ratio, taper, and sweepback are given in the appendix. 

Torsion. — The previous analysis has ignored the effects of 
wing twist due to torsion. The solutions obtained are, in 
reality, those for wings of infinite torsional stiffness. In 
general, since the flexural axis (or torsion center) is behind 
the center of pressure at all spanwise stations of the wing, 
the twist of the wing due to torsion will tend to compensate 
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for the twist due to bending. For wings having large angles 
of sweepbacli such as are necessary for efficient supersonic 
flight, the aerodynamic twist due to torsion has been calcu- 
lated to be about 15 to 20 percent of the twist due (o bending 
(for thin wings). In such cases, the effect of the torsional 
deformation on the spanwise loading may be neglected in 
calculating the torsional moment. Equation (7) may be 
utilized in the calculation of the torsional moment in this 
instance. A complex simultaneous solution is thereby 
avoided. 

An expression for the torsional moment at the root section 
of the wing beam (perpendicular to the elastic axis) may be 
obtained by assuming that the distance from the center of 
pressure to the flexural axis for any section of the wing is a 
constant percentage of the local chord . Then 

T l 

fi Y =S f (C0S A) Jo a 0 dy 

where c is the local streamwise chord given by the equation 
C=S Co j^l (1 X) 

and X denotes the taper ratio of the wing and f the distance 
from the center of pressure of the section lift to the flexural 
axis in terms of the streamwise chord. The function describ- 
ing the spanwise loading Ijq is given by equation (7). 

The equation for the torsional moment at the root may be 
written as 

(3 |- r == s 2 T c 0 cos A JJ ^ [i-d-x) |J dy 

p L = s * f c 0 cos A [JJ ~ dy-± (1-X)/; l - y dy ] 

As will be shown later, it is convenient to derive the ratio of 
the torsional moment at the root to the bending moment at 
the root. The bending moment at the root is given as 


and 


where 


and corresponds to the spanwise center of pressure for the 
load on the half wing. The value of Y may be determined 
by a mechanical or analytical intergration of equation (7). 

When the assumption is made that the twist due to torsion 
varies linearly across the span (or that the ratio T/6J is 
constant across the span), the incremental angle of attack 
of any section of the wing due to torsional deflection may be 
written as 

cos A Ts' Ts y 
aT ~ p GJ y ~lDp 



_Ms y ( T r l r E\ 
aT ~ El p \M T J t GJ 

and by adding this expression to the angle of twist due to 
bending (equation (1)) the total angle of twist of any section 

a ‘ =a ~m | [ tan A_ (tS) SI (14a) 

or 

a ‘ =a_2 i~ Wf3 [ tan A ~ S] (14b) 

Combined bending and torsion. — Expressions for the 
aerodynamic properties of swept wings experiencing both 
bending and torsional deformation may be obtained from 
equations (6) to (12) if tan A is replaced by 


The equation for the angle of attack at maximum load 
factor for combined bending and torsion is then 


0 2 2 ( S p .32 m 5 / 2 

~4m 2 Fi(2s 2 ^ + 3/3 3 tt (m+1) 

: m . )&]«} 


Pyiaz £ 

2 ^~d 




tan A— 


(15) 


where 


2V 

M T ' 


: fc 0 




cos 2 A 


In applying the foregoing analysis to a specific wing, it is 
convenient to use the equations to obtain the ratio of C„ a 
or C La for the elastic wing to the value for the rigid •wing. 
Multiplying this ratio by the value of C La or for the rigid 
wing as determined by the complete theory wherein the 
region within the Mach cone of the tip, and so forth, is 
considered, will then give more accurate parameters for the 
elastic wing. Then 


. I 4 m 2 s 
qa~ • 


a rigid 


a elastic ^ 


C m 


a rioid 


3j8 z 7t (m+l) 2 E d r 

8 2 Vw 

Gmaz £ 

3 jSx(m+l) 2 

E d,a n L 

8 2 Vto 

&max g f 

3 /Stt(to+1) 2 

E d T a n [ 


(16) 


(17) 


(I r E\ 7V-]E< 

\Z&) Mr] F s 

(18) 


In using the preceding equations, it is necessary to solve for 
a n . This, in turn, involves finding the ratio T T /M r which is 
determined by the parameter Y (usually has a value of about 
0.40). 

A solution of the combined bending and torsional deforma- 
tion effects can be obtained by assuming a value of Y, solving 
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for a n , and checking the value of Y from a moment and area 
integration of a plot of equation (16) to see if a second ap- 
proximation is required to determi n e a n more accurately. 

The previous equations apply primarily to fiat lifting wings 
or to twisted and cambered wings for which the loading due 
to twist and camber is the same essentially as the loading 
due to change in angle of attack. For wings with somewhat 
arbitrary camber and/or twist, these equations apply to all 
accelerated flight conditions. A solution for the aeroelastic 
characteristics in steady level flight for such wings must 
involve a consideration of the effects of the loading due to 
the known arbitrary twist. 

WING WITH A SUPERSONIC LEADING EDGE 


The foregoing analysis has treated wings with the leading 
edge swept behind the Mach cone. The same method, 
however, may be applied to wings swept ahead of the Mach 
cone. In this case, however, the expression for the pressure 
field for the incremental twist at any spanwise station, cor- 
responding to equation (2), is given by reference 10 as the 
real part of 


A p 4 a m 1—mt 

2 iSs-ym 2 — 1 \t—m\ 


(19) 


where a, t, and m are as defined for equation (2) . 

Expressions for the pressure distribution, lift, moment, 
and load distribution may be obtained in the same maimer 
as for a wing with a subsonic leading edge although the inte- 
grations a,re more involved. 

DISCUSSION 

SUPERSONIC LIFTING-SURFACE THEORY 

The results of the foregoing analysis are best illustrated 
by applying them to a specific wing. For this purpose, the 
wing shown in figure 4 was selected, having the geometric 



Figure 4. — Sketch of the wing used, in the calculations. 
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(a) Span load distributions for the rigid and elastic wings at the same angle of attack. 
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(b> Span load distributions for the rigid and elastic wings of the same total lift cofficienfc. 
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Figure 5. — A comparison of the span load distributions for the rigid and the elastic wings 
in accelerated flight. Jlf, 1.414; a-mu, 30,000 pounds per square inch; nW/S t 150 pounds per 
square foot; q, 211.2 pounds per square foot. 


and structural material characteristics given in the table in 
the figure. The calculations were made for various values 
of the parameter nW/S and for two values of the maximum 
design stress.* 

Span load distributions for the wing are shown in figure 5 
for a Mach number of 1.414, a value of vWjS of 150 pounds 
per square foot, a design stress of 30,000 pounds per square 
inch, and a dynamic pressure of 211 pounds per square foot 
which corresponds to flight at 60,000 feet altitude. The load 
distribution curves of part (a) of figure 5 are for the same 
angle of attack of the root section and show that the elasticity 
of the wing results in an appreciable decrease in lift-curve 
slope. In this case, the reduction experienced by the elastic 
wing amounts to 15 percent of the value for the rigid wing 
of the same plan form. Part (b) of figure 5 shows the load 
distribution curves for constant total lift coefficient. These 
load-distributions are of si gnifi cance in illustrating how the 
change in span load distribution due to elasticity may be 
expected to shift the longitudinal center of pressure forward. 

4 Calculations show that the wing has sufficient depth to withstand the maximum loading- 
assumed without failure. 
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COMPARISON OF AEROELASTIC EFFECTS AT SUPERSONIC SPEED WITH 
INCOMPRESSIBLE FLOW SOLUTIONS 


In. calculating lift and stability characteristics of elastic 
wings, it should be noted that errors resulting from assuming 
the extent of the wing beam as given in figure 1 and from 
ignoring the lift within the tip Mach cone may be minimized 
by using the analytical expressions which give the ratio of 
lift-curve slope or the ratio of moment-curve slope for the 
elastic wing to that for the rigid wing. These ratios may be 
used with the rigorous values of C ma and C La from reference 5 
to obtain accurate values of C m<t and C La for the elastic wing. 

Such ratios have been computed for the wing shown in 
figure 4 as functions of the dynamic pressure at a flight 
Mach number of 1.414. For comparison, the same ratios 
have been computed as functions of the dynamic pressure 
for incompressible flow by the theory of reference 11. Fig- 
ures 6 and 7 show the results of these calculations which 
were made for two values of_n W/S of 150 and 300 pounds 
per square foot and two values of design stress, 30,000 and 
45,000 pounds per square inch. Figure 8 shows the shift in 
neutral point 5 due to wing elasticity as calculated from the 
data of figures 6 and 7. 

The results indicate that the adverse effects of the aero- 
elastic deformation of the wing are a little more severe at 
supersonic speed. At constant dynamic pressure, the differ- 
ences in the aeroelastic effects as computed by incompressible 
flow theory and by supersonic lifting-surface theory are found 
to be due largely to the fact that the center of pressure of 
the sectional lift is farther forward at subsonic speed, result- 
ing in a difference in torsional deformation which compen- 

1 Neutral point is defined as the position of the center of gravity along the mean aerodynamic 
chord for neutral stability. 
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Figure 6. — Variation with dynamic pressure of the ratio of lift-curve slopes for the elastic 
and rigid wings in accelerated flight. 
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sates somewhat for the bending deformation. The compar- 
ison indicates that the dynamic pressure is the primary 
variable involved in determining the aeroelastic character- 
istics, at least for wings swept behind the Mach lines. 

In regard to the range of application of the equations, 
calculations made using more rigorous structural theory with 
simple strip theory show that the method of the present 
report may be expected to give accurate estimates of aero- 
elastic effects as great as, for instance, a 30-percent loss in 


lift-curve slope. "Within such limits it is expected that the 
estimate of the neutral point shift due to elasticity willj.be 
much more accurate than for analyses using elementary 
aerodynamic loading. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec. 3, 1943. 


APPENDIX 

mathematical derivation of loading functions and plan-form constants 


The functions /i (y) , j*(y) , f 3 (y) , and / 4 (y) and the constants 
F u F 2 , F 3 , Fi which appear in equations (7) to (18) of the 
text are given in this appendix. Analytical expressions for 
the constants F 2 , F s , F^ are found to be very long and tedious 
to use. It may prove easier to evaluate these constants by 
graphical integration of the corresponding integral equations. 

The functions fi(y) and F(y) were developed from the 
following integral: 
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The constants F t and F a are evaluated as follows: 
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Then 
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